Effective transition operators for resonant inelastic X-ray scattering 
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Effective symmetry-based transition operators for resonant inelastic X-ray scattering (RIXS) are 
derived that show how the scattering between different states depends on the polarization of the 
incoming and outgoing X-rays. In spherical symmetry, the effective operators can be rewritten in 
terms of spin operators, although the expressions depend on the nature of the ground state. For lower 
symmetries, the combined action of the crystal field and the spin-orbit interaction breaks up the spin 
space and spin operators are no longer appropiate operators. By taking iridium compounds as an 
example, it is demonstrated that effective scattering operators can still be obtained. These effective 
transition operators facilitate our understanding how RIXS couples to elementary excitations. 
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Collective excitations play an important role in the 
electronic and magnetic properties of materials. A typi- 
cal example is a magnon, a spin-flip excitation that dis- 
perses through the solid via the exchange interactions. 
These excitations are typically studied with inelastic neu- 
tron scattering where magnons are created through the 
interaction of the spin of the neutron with the magnetic 
moments of the material. However, magnon excitations 
can also be studied via the inelastic scattering of X-rays 
[ll-Q. The creation of electronic excitations can be en- 
hanced by tuning the X-ray energy to a particular res- 
onance where an electron from a deep-lying core level 
is excited into the valence states. Using this technique, 
known as resonant inelastic X-ray scattering (RIXS) Q, 
magnon dispersions throughout the Brillouin zone have 
been mapped for copper-oxide based materials [1, 
More recent experiments on iridates 0, § provide an ex- 
cellent example of how RIXS and inelastic neutron scat- 
tering can complement each other. Iridates have come 
to the forefront recently [ol-fl^ as a prototypical exam- 
ple of a system where the orbital and spin of the valence 
electrons are entangled by the strong spin-orbit inter- 
action. However, the small size of iridium-oxide single- 
crystals and its large atomic weight significantly hamper 
the study of these materials using neutrons. Using RIXS, 
Kim et al. @] were able to map a magnon dispersion with 
a width of about 0.2 eV in Sr2lr04. In addition, excita- 
tions were observed 0, Q with an energy of 0.5-0.8 eV 
which were assigned to entangled spin-orbital waves. 

One of the attractive features of inelastic neutron scat- 
tering is that its cross section is due to a comparatively 
simple dipole interaction between the neutrons and the 
magnetic moments in the solid. The RIXS spectral in- 
tensity, on the other hand, is due to a complex process, 
where after absorption of the incoming X-ray, the sys- 
tem is in a complicated intermediate state dominated by 
strong interactions with the core hole. However, it is 
exactly the large core-hole spin-orbit coupling and the 
many-body effects between the core hole and the valence 
shell that make spin flips possible. Subsequently, the flnal 



states are reached via the radiative decay of the core hole. 
Analytical evaluation of the resonant scattering ampli- 
tude [T3)-fl6| can provide more insight into the spectal 
weights of the flnal states, but generally require approx- 
imations of the intermediate-state propagator. However, 
when studying the momentum and polarization depen- 
dence of magnons and orbitons, one is less interested in 
the intensity compared to features of different symmetry, 
but more in the relative strength of the components of 
a particular symmetry. For example, for magnon excita- 
tions, it is important to understand how the polarization 
conditions affect the ratio between spin conserving and 
spin-flip transitions. In this Letter, it is shown that these 
transitions can be expressed in relatively simple opera- 
tors. Although the operators can often be expressed in 
terms of spin operators [13 , this result is not general and 
depends on the nature of the ground state. The relative 
strength of the operators can be tuned by the polariza- 
tion vectors of the incoming and outgoing X-rays. 

The RIXS intensity is given by the Kramers- 
Heisenberg equation 6], 
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where the energy loss h{u! — uj') of the photon equals 
the difference Ef — Ei between final and initial energies 
of the system; e and e' are the polarization vectors of 
the incoming and outgoing X-ray photons, respectively. 
The scattering amplitude consists of dipole operators D 
describing the absorption and radiative decay and an 
intermediate-state propagator: 
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where H is the Hamiltonian of the system and F the 
intermediate-state lifetime broadening. In order to treat 
systems with different local symmetries, we define |m) — 
\^mlm) for m = i, f , where F^ (not to be confused with 
the lifetime broadening) is the total symmetry and 7m its 
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component. For example, for a divalent nickel ion in oc- 
tahedral symmetry (Oh), the total symmetry is F —'^ A2 
and the components 7 — 1, 0, —1 are the three spin states 
of the triplet. After recoupling , the scattering am- 

plitude, 
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k=0 q 

can be separated into a simple numerical factor au = 
{2k + l)nl^^ [11], a term Tf{e,e') = [e'*,e]J describing 
the dependence on the polarization vectors as a tensor 
product [1^ , where q = F7 is the symmetry of the tran- 
sition operator, and an effective scattering operator F^. 
Using the Wigner-Eckart theorem [23|, we split the ma- 
trix element of i^^ into two terms 

(/IF^^'MK) = {Tf\\F'^{u;)\\r,)n,MfW), (3) 

where the reduced matrix element (F/||F'^(a;)||Fi) de- 
scribes the strength of the transition between the man- 
ifolds of symmetry F^ and F / . The second term can be 
written as an effective operator 

W, = J2\f){m){^n-^„ (4) 

where fk are factors that are introduced to obtain more 
nicely defined operators. The term {f\qi) is a coefficient 
that gives the coupling between the initial state i and the 
final state / via a tensor with symmetry q. It is impor- 
tant to realize that all the dependence on the interme- 
diate states is contained in the reduced matrix elements 
that only depend on the total symmetries F^ and F^. The 
polarization dependence of the intensities of the different 
components of a particular manifold only depends on the 
effective operators Wg. Although the above result is ex- 
act, it is not directly clear that the operators Wq are of 
any practical use since they depend on coupling constants 
{f\qi) which are matrix elements between complex initial 
and final states. In the following, it is demonstrated that 
the Wq can be expressed as relatively simple operators. 

Since the nature of the operators Wq depends on the 
symmetry of the problem, we focus on some specific situ- 
ations. Let us first consider spherical symmetry. Typical 
systems are rare-earth ions, where generally the orbital 
moment L and spin S are coupled to a total moment J 
via the spin-orbit coupling. The total symmetry of the 
state is characterized by the total momentum (F,„ — >■ J); 
the components are given by the projection (7™ — Mj). 
This gives \i) = \JM,j) and |/) = \J'M'j). In a mag- 
netic field along the positive z axis, the ground state has 
Mj = — J. For the scattering operators, all the compo- 
nents of the scattering term of rank k belong to the same 
representation and F7 — >■ kq with q = fc, fc — 1, • • • , —k. 
The coupling constants in Eq. Q are directly related 
^ to the Clebsch-Gordan coefficients {J' M'j\kqJMj). 
Obviously, well-defined effective operators are the compo- 
nents j'^p = Jq of the total angular momentum J, whose 



matrix elements are proportional to {J'Mj\lqJMj), and 
combinations thereof, such as JqJq', that give rise to 

(k) 

higher-order terms Jq . Expressing the operators Wkq in 
terms of spin requires some more care, since the ground 
state is not directly an eigenstate of the spin. Let us con- 
sider excitations within the ground-state manifold char- 
acterized by J. If the ion is coupled via exchange inter- 
actions to other magnetic ions in the lattice, these exci- 
tations would lead to spin waves. For these states, we 
can define general spin operators of rank k via 

= J2 \SS',){SS'JkqSS,)nslk{SS..\, (5) 

where S is the spin. <t*^°) is the identity operator; the 
Pauli spin matrix is given by (t^^^=S/S; the spin op- 
erator of rank 2 is normalized by cr'-'^>=S'''^y S{2S—1). 
Note that the maximum change in spin is given by k. In 
spherical symmetry, we can relate the spin operator to 
the total angular momentum operators Jq via 

(JM'jlcT^J'^lJMj) = ALSJk{JM'jl4''^\JMj), (6) 

where the A^sjk gives the details of the coupling between 
the orbital moment L and the spin S in the ground state 

ALSjk = V(2J+1)(25+1)^ [s S l]' 

where the bracketed term is a 6j symbol. When J = 
L + S, one has A^sjz = 1 and the "spin-only" result of 
Haverkort [l7| is reproduced. However, when the spin 
is antiparallel to the orbital moment, the weighting of 
the different (tC') to the resonant X-ray scattering ampli- 
tude changes. For example, for J = L — S* with L > S, 
Alsji = -J/(J+l)and^isj2 = J(2 J- !)/( J+l)(2 J+ 
3). Therefore, in spherical symmetry, the inelastic scat- 
tering amplitude can be expressed in terms of spin oper- 
ators, since the spin operator is a proper spherical tensor 
However the expressions depend on the coupling between 
the orbital and spin moment in the ground state. 

An interesting feature of scattering within the ground- 
state manifold is that the reduced matrix element can be 
expressed in terms of the elastic scattering amplitude 

J-iH = {J,^J\Fl;{uU-J) = {J\\F''{uj)\\J)njjk. 

which is directly related to the X-ray absorption, which 
is proportional to — Im[J"gj (w)]. For the ground-state 
manifold with total angular momentum J, the inelastic 
scattering from Mj to M'j, of importance for magnetic 
excitations can be written as 

Tico, e, e') = ^ ^Ti{u.)[e'*,ef ■ (8) 

where the polarization dependence is directly coupled to 
the effective operator. 

The correspondence between the effective operators 
and the local spin operators holds in many cases, but 
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fails for lower symmetries where the combined action of 
the crystal-field and the spin-orbit coupling can effec- 
tively break up the spin space. A representative example 
is a low-spin tetravalent iridium ion. Iridates have re- 
ceived significant attention recently i7h121| since in these 
systems, despite a large bandwidth, the spin-orbit inter- 
action is not quenched. This can lead to new kinds of 
topological orderings. has a configuration (^T2 

in Oh)- Orbital and spin space are coupled via the spin- 
orbit interaction and the symmetry splits into E" and 
U' symmetries, which are two- and four-fold degenerate, 
respectively. For an iridium ion, the E" representation is 
the lowest in energy. This still holds in planar symmetry 
{D4h) and the wavefunction for the t2g hole in the ground 
state can be written as 

|£;",±i) =sina|0,±i) -cosa|±l,Ti), (9) 

where |0) — —i\xy) and | ± 1) = -^{—\zx) ± i\yz)). Fig- 
ure [DJa) shows the x-ray absorption and x-ray magnetic 
circular dichroism at the L3 edge. The spectra are cal- 
culated using standard multiplet codes including a 
cubic crystal field of 3 eV. The symmetry is further low- 
ered by applying a planar crystal field, which also de- 
termines the relative occupation of the xy and yz/zx 
orbitals. The spectra are relatively simple: the main line 
is due to the four empty Cg states; at the low-energy 
side, there is intensity due to the empty t2g state, which 
shows a dichroic signal. Let us consider the scattering 
for fc = 1 between the different Sz components of the 
ground-state manifold \E"Sz) with 5'^ = ±i. In planar 
symmetry, the three components of the k = I transi- 
tion operator are r7 = A2; E,l; E, ~1 which are given 

by Wa, = -jai^^ = -^Sz and We,±i = ^^f ^ = jS,/y 
where the spin operators are given in Eq. ([5]) using the 
equivalence between E" and S = ^. Although the oper- 
ators are relatively simple, it is important to note that 
they are effective spin operators (indicated by the over- 
line in Si) between the two components of the E" ground 
state and not real spin operators. This becomes obvious 
when considering the limit a — > in Eq. © where the 
low-energy manifold is given by \E",zt^) — \zLl,^^). 
The matrix elements of the spin operator are given by 
{E"S'z\S^\E"Sz) = -Sz6s^,s'Az. Therefore, from the 
point of view of the (real) spin, the \E",±^) manifold is 
an effective Ising spin due to the fact that a spin fiip must 
be accompanied by a change in orbital. This implies that 
if the inelastic x-ray scattering amplitude would be de- 
scribed solely in terms of spin operators, magnon excita- 
tions between \E"^) and \E",—^) would not be allowed 
which clearly contradicts the finite spin-fiip scattering 
found numerically, see Fig. [IJa). 

For excitations within the manifold with the lowest 
energy, the fc = 1 transition operator is proportional to 

J2 iE"\\Fh{u;)\\E")pM'^ - {E"\\Fliu;)\\E")Pzai'^ 
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FIG. 1: (a) The imaginary part of the X-ray scattering 
amplitudes [-{l/-K)lm.{E" s^\Fr^\E" \)] for an Ir^+ ion {ti,g). 
Shown are the isotropic X-ray absorption (A: = 0, = i), 
the x-ray magnetic circular dichroism (XMCD) (fe = 1, = 
5), and the inelastic spin-flip scattering amplitude (fc = 1, 
= — I) in planar symmetry {D^h) in the limit that the 
hole density in the xy orbital n^^ goes to zero. For compari- 
son, the inelastic spin-flip scattering in octahedral symmetry 
is shown, (b) The reduced matrix element for the k = 1 tran- 
sitions as a function of the hole density n^^ in the xy orbital. 
The terms with and without spin flip are given hy T — E, A2, 
respectively. There are no spin flip terms at the L2 edge. 

with Pi = Tl{e,e') = (e'* x e)i with i = x,y,z. Within 
planar symmetry, the Sz and spin-flip terms have a dif- 
ferent symmetry and their relative weight therefore de- 
pends on the frequency-dependent reduced matrix ele- 
ment. Since the spin excitations occur in the ground- 
state manifold, we can express the inelastic scattering 
amplitudes in terms of the elastic one. For the A2 rep- 
resentation, (£;"||F;^^(a;)||£;") = FX^_^i {uj)/nii^, where 

— Im[i^^^ ^[ (w)] is the XMCD with the magnetic or ex- 
change fields in the z direction. However, the cou- 
pling coefficient is zero for E symmetry in this config- 
uration. It is nonzero when the magnetization is in 
the plane, {E"\\FX^{u;)\\E") = Fi;,_,,(^)/'^iii where 
Im[Fj;^ (cj)] is the XMCD measured along the x di- 
rection. The values of the reduced matrix elements are 
shown in Fig. [ijb) as a function of the number of holes 
n^y in the xy orbital. Although there are no spin-flip 
terms at the L2 edge, the spin-conserving and spin-flip 
terms are of comparable order of magnitude at the L3 
edge. When there is no planar distortion, (n^^ = 0.35), 
all directions become equivalent and the scattering am- 
plitude becomes proportional to P ■w'^^\ Since there are 
no spin flips at the L2 edge, this also implies that the 
reduced matrix element for ai^"* must be zero, see Fig. 

The approach described here is not restricted to ex- 
citations within the ground-state manifold, but can also 
describe transitions between different manifolds. Let us 
consider again the situation for an ion and take oc- 
tahedral symmetry [cos^ = | in Eq. ^] . The wave- 
functions for U' symmetry are |/7',±|) = =F|±l,±i) 

and |J7',±i) =±y||0,±i)±yi|±l,=F5). Let us con- 
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sider the case of fc = 1, which in octahedral symmetry 
corresponds to the Ti representation with components 
q — 1,0, —1. The effective transition operator given by 

Wr.g - \U\S,+q){U',S,+q\T^q,E"S,){E"S,\ (10) 



brings the system from the ground state with Sz 
into the four-fold degenerate final states. The matrix 
elements of Wti? can be directly related to a Clebsch- 
Gordan coefficient 



{U',Sz+q\T^q,E"S, 



, Sz+q\lq, \Sz 
2m 1 3 




and the operator can be written as an effective total 
angular momentum J. The contributions of Wtkj to 
the scattering amplitude can be varied via the polariza- 
tion dependence giving a weight —{Px + iPy)/^/2, {P^ + 
iPy) I \f2,Pz for q = 1,0,-1, respectively. This allows 
one to probe different components of the orbiton state. 



In conclusion, it has been shown that when considering 
the polarization dependence of excitation of a particular 
symmetry, such as magnon or phonons, relatively simple 
effective transition operators can be derived that depend 
only on the local symmetry. The coupling of these op- 
erators to the polarization dependence allow the tuning 
of excitations into the different components of the final 
state (for example, transitions into states with and with- 
out spin flip). 
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